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Abstract 

We discuss implications of Planck results for models with local type non-Gaussianity. 
In light of the recent results of the Planck satellite, we constrain model parameters 
of several representative models and give the prediction of trispectrum, in particu- 
lar, gNL- We also consider interesting possibilities that trispectrum appears as the 
first signature of the non-Gaussianities of the curvature perturbations, that is, /nl 
is small while (/nl can be significantly large. 



1 Introduction 



Very recently, the Planck mission has released data from the first 15.5 months of Planck 
operations for cosmic microwave background (CMB) anisotropics pQ. They determined 
cosmological parameters with unprecedented accuracy such as the baryon, the dark mat- 
ter, and the dark energy densities, which strongly support the so-called concordance model 
of cosmology [2]. They also gave strong constraints on primordial curvature perturbations. 
The spectral index n s = 0.9603 ± 0.0073 (68%CL) [3] significantly deviates from unity, 
supporting the slow-roll inflation paradigm. Any deviations from Gaussianities of primor- 
dial curvature perturbations are not found. In particular, the local type of /nl is now 
strongly constrained as —8.9 < /nl < 14.3 at two sigma level [I], which rules out a lot of 
light field models predicting large local type non-Gaussianities. 

In fact, one may wonder if light field models such as the curvaton [HH7] and the mod- 
ulated reheating [HI [9] scenarios might be excluded because they are often claimed to 
generate large non-Gaussianities. However, this is not the case. For example, as is well 
known, in the curvaton scenario with a quadratic potential, the local type of /nl is given 

by USE] 

5 5 5r 

4r 3 6~ (1) 
Here r is roughly the fraction of the curvaton energy density at the curvaton decay and is 
defined as 

3p CT 



(2) 

decay 



3p a + 4p. 

where p a is the curvaton energy density and p r is the radiation energy density. The 
constraint -8.9 < / NL < 14.3 (95% CL) is equivalent to 0.08 < r (< 1) (95% CL), which 



rules out the curvaton model with small r F 1 ! . However, to be fair, a natural value of r 



without fine-tuning is unity because the curvaton easily dominates the energy density of 
the Universe since the curvaton behaves like matter while the other components behave as 
radiation. Such a value of r(= 1) yields /nl = —5/4, which is still allowed by the recent 
Planck data. Thus, a simple and natural model of the curvaton is still viable. 
In the same way, the modulated reheating scenario predicts /nl as 



/nl = 5 (l - , (3) 

where V is the inflaton decay rate depending on the modulus cr, T a = dT/da, T aa = 
d 2 T/da 2 , and the inflaton is assumed to oscillate around its minimum with a quadratic 
potential. The constraint -8.9 < / NL < 14.3 (95% CL) leads to -1.9 < TT aa /Tl < 2.8 
(95% CL), which suggests that the second derivative of T with respect to cr is significantly 



^Here, we simply adopt the Planck constraint on /nl to the constraint on the model parameter r. 
However, a direct likelihood analysis of r with adopting a prior < r < 1 performed by Planck team gives 
the constraint 0.15 < r (95% CL) [4]. This is the reason why the constraint by the Planck team is a bit 
different from ours. We thank Nicola Bartolo for the clarification on this point. 
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suppressed. Thus, as long as T linearly depends on the modulus a, /nl is predicted to be 
5, which is still allowed by the recent Planck data. 

Thus the light field models are still viable. After excluding the model parameters 
inconsistent with Planck data, we end up with the light field models that generically 
yield /nl = Given that the standard inflation model, in which the inflaton itself is 

responsible for the curvature perturbations, predicts (9(0.01) [HUH)], it is still important. 
Then, it is still important to detect order of unity /nl for discriminating light field models 
from the standard inflation model though, in this case, we have to take into account the 
intrinsic CMB bispectrum /nl < 1 coming from the second order effects of the evolution 
of the curvature perturbations [T6lfT7]. 

In this paper, we first discuss how much light field models are constrained according to 
the Planck data. As explained above, both of the curvaton and the modulated scenarios 
are still allowed, and, in some sense, the constrained value of /nl is reasonable to avoid 
fine-tuning. Then, we give constraints on model parameters of light field models including 
the curvaton and the modulated scenarios and predict the trispectrum, in particular g N L, 
for each model. 

Next, we are going to pursue another interesting possibility. Though the recent data 
of the Planck satellite claims that bispectrum of the curvature perturbations is small, it 
does not necessarily imply that the non-Gaussianities of the curvature perturbations are 
insignificant because such non-Gaussianities might first appear on their trispectrum. In 
the latter half of this paper, such a possibility will be discussed in detail. Actually, for 
example, if we consider the modulated reheating scenario and its modulus also behaves 
like a curvaton, /nl can be almost canceled and be small, but #nl still can be large. We are 
going to discuss such possibilities and investigate which combination of model parameters 
can realize such possibilities. 

The organization of the paper is as follows. In the next section, we constrain the 
model parameters of each light field model based on the recent Planck results and give 
the predictions of the trispectrum $nl, whose observations are essentially important for 
pinning down the model. In Sec III. we discuss models, in which /nl is small, actually, 
within the constraints given by the Planck satellite, but $nl can be large as the first signal 
of the non-Gaussianities of the curvature perturbations. Detailed explanations of why <?nl 
can be large while keeping /nl small for such models are given. The final section is devoted 
to summary of this paper. 

2 Light field models 

During inflation, there can be light fields other than the inflaton. They acquire quan- 
tum fluctuations during inflation, which can be converted to the curvature perturbations. 
Though a lot of conversion mechanisms have been proposed [18], the S N formalism [T9H23] 
enables us to make a systematic treatment to evaluate the final curvature perturbations. 
According to the SN formalism, the superhorizon curvature perturbations at the final time 
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t = tf can be easily estimated by 

C(t f ) = N a S V : + l -N ab 5 V a M + , (4) 

where t* is the time shortly after the horizon exit, a subscript a, b, and c represents a 
light field, N a = dN/dip®, and so on. Here, 8tp% represents a field fluctuation evaluated at 
t — t* and is assumed to be Gaussian. Then, the bispectrum and the trispectrum of the 
curvature perturbations are characterized only by the three parameters, /nl 3,1 , g^fi an d 
TN°£ al as follows, 

(CfeCfeCfe) = (27r) s fl c (fci, fe, + + h), 

(CkSk^kJ = ( 2 ^f T d k ^ h, h, h)S(h + h + k 3 + k 4 ), (5) 

where 

B c (ki, h, h) = ^/N L cal (Pc(h)P c (k 2 ) + P c (h)P c (h) + P c (ks)Pdh)) , (6) 

?c(A;i, k 2 , h, h) = r^ 1 (P c (k 13 )P c (k 3 )P c (h) + 11 perms.) 
^4 

+ -^ al (Pc(k 2 )P C (h)P c (h) + 3 perms.) , (7) 

with &i3 = | fc x + fc 3 1 . These three parameters are easily evaluated at the tree level according 
to the 5N formalism, 



6 local N a N b N ab 

5 /NL " (Ayv<f' 

local iV a iVyV ac iV c b 



(N d N d Y 



54 _ AWiV a iV fe iV c 

25 " (N d N*) 3 • (8) 

Below we omit the suffix "local" for simplicity. There is no general relation between /nl 
and #nl because (/nl depends on the third derivative of N in addition to the first one. On 
the other hand, there is a general inequality between /nl and tnl [21], 

TNL > (lfw) , (9) 

When only one source (field) contributes to the curvature perturbations, the equality must 
hold. On the other hand, when there are multiple sources, the equality can be violated. In 
the following analysis, we concentrate on single source case for simplicity. In this case, tnl 
is completely determined by /nl as tnl = 36f^ L /25. Therefore, #nl is a key observable 
quantity to discriminate light field models. Below, we constrain model parameters of light 
field models and give the prediction for g^h based on the recent Planck satellite results. 
The extension of our discussions to multiple source case is straightforward and, generally 
speaking, the allowed region of the model parameters is widened. 
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2.1 Cur vat on 



In the curvaton model, the curvaton fluctuations are converted into the curvature pertur- 
bations when the curvaton decays into relativistic degrees of freedom, which occurs after 
inflation. The important quantities determining the resultant curvature perturbation are r 
representing the curvaton fraction to the total energy density and tr osc , the curvaton value 
when the curvaton starts oscillations. If a curvaton potential deviates from a quadratic 
form, (j osc generally depends on the curvaton value cr* at the time of horizon crossing. The 
non-linearity parameters from the curvaton are given by [25] 



/nl4 1 + ^= "2-r, (10) 



6 „ 3 /„ cr osc a^ 

5— 2r V" ' < 2 SC 

1 I . „ CoscCqsc \ 



2 V a. 



12 

osc 



+ 10r + 3r z , 11 



where tr osc = do- osc /da* etc. Though the relation between a osc and tr* is nontrivial in 
general, <7 osc has a linear dependence on a osc for a quadratic potential of the curvaton. In 
this case (<r ' sc = <7 g C = 0), the non-linear parameters reduce to 

§An. = | - 2 - r, (12) 
54 Q 1 

^NL = -- + - + 10r + 3r 2 , (13) 
which lead to the following consistency relation, 



^"54 



54/ 2 L - 60/ NL - 125 - (9/nl + 5) a/36/ 2 l + 120/ NL + 250 



For /nl ^1 (r < 1), this relation is approximated as [18 



(14) 



1 575 

Snl = - Y /nl- — + 0(r). (15) 

Thus, the present constraint on / NL , —8.9 < /nl < 14.3 (95% CL), leads to 0.08 < r < 1 
(95% CL). The allowed region of #nl is predicted to be —53 < #nl < 2 (95% CL). 



2.2 Modulated reheating 

In the modulated reheating model, the decay rate of the inflaton T depends on some 
light field called modulus o . After inflation, the inflaton starts the oscillation around its 
minimum. When the potential around the minimum is well approximated by a quadratic 
type, the energy density of the inflaton oscillation decays in proportional to a~ 3 (a: the 
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scale factor) and hence the oscillation behaves like a non-relativistic matter. On the other 
hand, the energy density of the Universe after reheating is dominated by radiation, whose 
energy density decays in proportional to a -4 . Thus, the fluctuations of the decay rate of the 
inflaton leads to the energy density (curvature) perturbations of the Universe. Assuming 
that the decay rate is much smaller than the Hubble parameter evaluated at the end of 



inflation |*3> the final curvature perturbation is easily evaluated by the 5N formalism as 



C = --—8(7* + — (-^ + + J_ f + 3^4^ - 2S ) 5al (16) 

s 6 r 12 V T V 2 ) * 36 V r 2 r 3 ' 

which yields the non-linear parameters, 

§/ N l = 6-6^, (17) 

54 / it r 2 r \ 

|^NL = 36 (2 - 3^ + . (18) 

It should be noticed that, when V linearly depends on a and the other higher derivatives 
vanish, /nl and #nl are predicted to be /nl = 5 and #nl = 100/3. 

Let us consider the following a dependence on T as a more concrete example |13j . 

a n 2 



a— + B—\, (19) 
M f M 2 ' 



where a and (3 are constants and M is some energy scale. The non-linear parameters are 
rewritten as 



Thus, the present constraint -8.9 < / NL < 14.3 (95% CL) leads to -0.9 < (3 /a 2 < 1.4 
(95% CL). In this example, the third derivative of T, T aaa , is negligible, which leads to 
the following consistency relation, 

50 

<?nl = 10/nl - y • (21) 
Therefore, ^nl is predicted to be —106 < #nl < 126 (95% CL). 



2.3 Inhomogeneous end of hybrid inflation 

In hybrid inflation, the inflationary phase are kept thanks to the positive effective mass 
squared m 2 of the waterfall field x and ends at the critical value cr of the inflaton <j) due 
to the tachyonic instabilities. Then, if the effective mass squared of the waterfall field 

# 2 This assumption does not necessarily hold in general. For instance, see [26] . 
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depends not only on the inflaton but also on another light field a, the critical value 4> cr 
also fluctuates, which leads to the perturbation of the duration of the inflation, that is, 
SN [27H31]. 

Following Ref. [29], let us consider a potential of the form, 



111 1 

i 1 2i2 2,- L 2 i2 , L r2 2 2 , L 22 



(22) 



which yields the effective mass squared of the waterfall field m 2 and the critical value of 
the inflaton <b cr as follows, 



2 , 2i2 , j}2 2 

-v +g <f> +f a , 



y/v 2 - /V 2 
9 



(23) 



Here /, g, and A are coupling constants and v is some scale related to the vacuum expec- 
tation value. The total duration (e-folding number) of the inflation is easily estimated 

as 



N 



1 



V_ 



(24) 



which fluctuates due to the perturbation of <p cr and generates the curvature perturbations, 



dN d4 a 1 
d(j) CT da * 2 



d 2 N 



da 



+ 



dN d 2 c 



da 2 



1 

+ 6 



d 3 N 



da 



d 2 N 



da 



d 2 



da 2 



dN d 3 c 



<90 cr da 3 



5ai 



Then, the non-linear parameters are given by 

6 
5 
54 



NL 



:9nl 



-M Pl V2e 
-(2e c 



cr XI 2 ' 
cr 

1 



/ NL + 2M^e c 



cr 



3 ' 



(25) 

(26) 
(27) 



25*"" v " 5 

where e = M^iV^/V) 2 /2 and 77 = M^V^/V are the standard slow-roll parameters and we 
have omitted the contributions comparable to the slow-roll suppressed parameters. The 
subscript "cr" represents the values evaluated at the critical point (p cr . From Eq. (T231 . we 
can easily evaluate the non-linear parameters as 



6 v 2 



54 2 t, 2 



(2? 



Then, the present constraint -8.9 < / NL < 14.3 (95% CL) leads to -10.7 < rj^v 2 /(fa 2 ) < 
17.2 (95% CL). We also have the following consistency relation, 



9nl 



10 t 



(29) 



Since r] cr < 1, #nl is predicted to be —30 < #nl < 48 (95% CL). 
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2.4 Inhomogeneous end of thermal inflation 

In the inhomogeneous end of thermal inflation model [32], the effective coupling between 
a flaton field and the cosmic temperature, g, depends on some light field called modulus 
a. During a mini-inflation phase, so-called thermal inflation, the flaton field is trapped 
at the false vacuum due to the thermal mass. When the temperature decreases down to 
the critical temperature depending on g, the flaton field starts to roll down to its VEV 
and the mini-inflation ends. Hence, the fluctuation of the effective coupling g leads to 
the fluctuations of cosmic e-folding number corresponding to the primordial curvature 
perturbations. The final curvature perturbation is given by 

c4-*-*+U— -4W+^(— -3^ +2 4W (30) 



2 g ' 4 V 9 9 2 J * 12 V 9 9 2 9 3 / 

where the subscript a denotes the derivative in terms of a. This yields the non-linearity 
parameters as 

^ NL = -2 + 2^p, (31) 

5 9a 



54 a n n99aa , 9*9 



—g NL = 42- 3—2- + — • 32 

25 V 9l 9% J 

In a similar way to the modulated reheating case, when g linearly depends on a and the 
other higher derivatives vanish, the non-linearity parameters are predicted to be /nl = 
—5/3 and #nl = 25/9. When the effective coupling, g, has the following a dependence as 

3 = 9 "{ 1 + a j} + ^)' (33) 

the non-linearity parameters are given by 



+ |— ,34) 



In this case, the present constraint —8.9 < /nl < 14.3 (95% CL) leads to —2.2 < (3/a 2 < 
4.8 (95% CL). 

In this model, we have also the consistency relation between the non-linearity param- 
eters as 

10 f 50 (k\ 

9~NL — — — JNL - [OO) 

where we have neglected g aaa and this leads —50 < gNL < 28 (95% CL). 
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2.5 Modulated trapping 

When the inflaton has a non-trivial coupling to other fields, the resonant particle produc- 
tion can happen. Such particle production significantly decreases the speed of the inflaton 
due to the backreaction effects. Then, if such particle production process depends on an- 
other light scalar field a through a coupling constant and/or a resonant point, the cosmic 
expansion is perturbed due to the perturbations of the light scalar field, which generates 
the curvature perturbations. 

Following Ref. [33], let us consider the coupling between an inflaton and fermionic 
fields x given by, 

£ int = --M(m - \(f))xx, (36) 

where M is the number of species of x particles with the same mass. Here, m and A are 
a coupling constant and the bare mass of Xi both of which are assumed to depend on 
another light field a. When the inflaton reaches the particle production point pp = m/A, 
the effective mass of x vanishes so that x particles are resonantly produced. The produced 
number density is estimated as 

\3/2 

™ P P = ^3l<U 3/2 > (37) 

where the subscript "pp" represents the quantities evaluated at the particle production 
time. Then, the equation of motion for the inflaton is modified as 

+ 3H ^ + ^M = MXn p J—) Q(t-t pp ). (38) 

u(p \ a pp J 

In order to quantify the particle production effects, we define 

A0(t) =0(t,A^O)-0(t,A = O), (39) 
which is easily evaluated from the equation of motion, 

A ^f A ^^f- (40) 

Here we regard H and dV (<fr) / deft almost constant because the duration of particle pro- 
duction is assumed to be short. According to the SN formalism, the final curvature 
perturbations are estimated as 



C = AN^ = -H^ = AV Jiy /2 , (41) 
0pp 187r 3 #p P 

which can expanded with respect to <5cr* as 

C = (AN(W) Ja* + \ (AN^) ac (5a*f + ± (AA^)^ (5a^. (42) 



S 



Though the general formulae for the non-linear parameters are a bit complicated, they are 
simply written in the case that both m and A are proportional to a, 




5e/3' 
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(43) 



54 



27 



(44) 




25e 2 / 9 2 ' 



where (3 is the "efficiency factor" given by 




(45) 



Then, the present constraint -8.9 < / NL < 14.3 (95% CL) leads to (3 < 0.10 (95% CL). 
We also have the following consistency relation, 



Then, g NL is predicted to be < #nl < 45 (95% CL). 
2.6 Velocity modulation 

In Ref. [31], it was shown that if particles a in the early universe have velocity fluctuation 
on large scales, their decay rate also acquires fluctuation through the fluctuation of the 
Lorentz factor and hence the curvature perturbation is generated via the same mechanism 
as the standard modulated reheating scenario. One of the scenario to realize such a 
situation is to assume that parent particle £ which decays into the daughter particle 
a has fluctuation of its mass <5ms due to its dependence on the light field having long 
wavelength fluctuations. Notice that because of the mass fluctuation 5ms, the resultant 
u particles receive not only velocity modulation but also density perturbation 5p a . As a 
result, the curvature perturbation generated at the time when a decays consists of two 
components, namely, the one generated by the standard curvaton mechanism and the 
other generated by the velocity modulation. Since both two components originate from 
the same fluctuations 5m-£, they are fully correlated to each other. Therefore, taking both 
the two effects mentioned above into account, the final curvature perturbation can be 
expanded in terms of 5m-£ as 



where each coefficient Ai(i = 1,2,3) depends on the four parameters in the model, Wq 
(equation of state parameter of a at the time of its generation), f2 s (density parameter of 
£ at the time of its decay), w a (equation of state parameter of a at the time of its decay) 



#NL = 




(46) 




(47) 
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and Q a (density parameter of a at the time of its decay). The expressions for the expansion 
coefficients A4 are lengthy (especially for A 3 ), we refer the readers to Ref. [31] for their 
explicit expressions. Assuming wo = ~ and f2s <C 1 just for simplicity, the non-linearity 
parameters are given by 



./j 



20 



NL 



^(36^(2^ - 1) + 5) + - 3w a tt a ) 



9nl 



9(3^(4^ - 1) + l) 2 Q CT ((3w CT - l)Q a + 4) _ 

+2(^(3^(24^(6^ - 7) + 35) + 10) - 7)fi CT + 2w a (9w a (8w a (6w a + 5) - 21) + 22) + 2% 

800 

2A3(3w CT (Aw a - 1) + iyttl(-3w a {l a + 4) 2 
(3w a - l)(w a (3w a (27w a (Aw a (9Qw a (2w a - 3) + 127) - 75) - 43) - 97) + 41)^ 

-3{w„(3w a {9w a {w a (l2w a {%w a (?>§w a + 19) - 413) + 2579) - 416) + 262) - 500) + 133)^ 
-2(^(36^(3^ - l)(4w a - 1) + 1) + l)(n a - 3wM„) 4 

-2(w a {3w a {3w a (3w a {12w a {12(37 - I2w a )w a + 133) - 5177) + 4570) - 728) + 674) - 517)fi CT 



+4( W(T + l)(9w a (w a (12w a (36w a (Aw CT + 9) - 181) + 265) - 2) + 209) 



We can easily find that /nl is at least 0(1) for any choice of (w a , f2 CT ). For fixed value of Q a , 
/nl takes a maximum ~ 15/f2 CT for w ~ 0.1. Numerically, we find that 0.5 < (7nl/tnl < 1 
for any (w a , Thus, a consistency relation between /nl and #nl can be written as 



n -r 36 f 2 
yNL — vm 25 ■' NL ' 



0.5 < C vm < 1. 



(50) 



Using the Planck result —8.9 < /nl < 14.3 (95% CL), the most conservative limit on <?nl 
in this model is given by < g NL < 294 (95% CL). 



(49) 



3 Trispectrum as first signature of non-Gaussianity 

In this section, we discuss interesting possibility in which trispectrum appears as a first 
signature of non-Gaussianities of the curvature perturbations. That is, we discuss the cases 
that /nl is small, as shown by the Planck results, but (?nl can be significantly large. Note 
again that tnl is also small for single source because we have the relation t nl = 36/ NL /25. 

Generally speaking, in order to realize large #nl while keeping /nl small (that is, 
N aa /N 2 = 0(1), N aaa /N% 3> 1), we can take two options. The first option is to consider 
large third derivative of N, that is, N aaa . Such examples include curvaton with a self cou- 
pling, modulated reheating model with the non-trivial (cubic) dependence of a decay rate V 
on the modulus a, and inhomogeneous end of thermal inflation with the cubic dependence 
of a coupling g on the modulus a. In these examples, as long as such cubic dependences 
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are extraordinarily large, #nl can be significantly large, which is quite manifest from the 
formulae given in the previous section. 



The second option is to make the first derivative No- accidentally smallj^j which leads 
to large #nl unless N aaa is also suppressed. It should be noticed that the second derivative 
N aa needs to be mildly suppressed to keep /nl of the order of unity. Such accidentally small 
N a can be obtained by tuning model parameters. For example, in the modulated decay of 
the curvaton model [35T437] . (essentially) two model parameters (r and YY^/Y"^) appear. 
By taking their adequate combination, N a (N aa ) is significantly (mildly) suppressed, which 
leads to large #nl with small /nl without resorting to N aaa . 

Such a situation also happens in the case that the same light field a contributes to the 
curvature perturbations multiple times in different ways. For example, in light field models 
discussed in the previous section, the fluctuations of a light field other than the inflaton are 
converted to the curvature perturbations. However, except the curvaton (and the velocity 
modulation), the evolution of such a light field is followed only until the conversion and 
the subsequent evolution is simply assumed to be negligible. However, generally speaking, 
such a modulus has only weak (gravitationally suppressed) interactions so that it is long- 
lived and can easily contribute to the energy density of the Universe at late times. That 
is, it is probable that a modulus in a light field model also behaves like a curvaton at late 
epoch, which becomes another origin of the curvature perturbations. In this setting, the 
curvature perturbations ( consists of two parts. The first part comes from each light field 
model contribution and the second part arises from the curvaton contribution, 

C = (light + Ccur- (51) 

It should be noticed that both contributions flight and Ccur arise only from the same pertur- 
bations Then, when the linear and quadratic terms cancel adequately but the cubic 
term does not, /nl becomes small while significantly large #nl can appear without resort 
to the cubic dependence of the couplings. This is exactly the second option we can take. 

In this section, we discuss this kind of possibilities in detail. More concretely, we 
investigate, for each model, which conditions on model parameters are necessary to realize 
such a possibility. 



3.1 Modulated Decay of the Curvaton 

Now in this section, we consider a model where the decay rate of the curvaton is modulated 
due to fluctuations of some other light scalar field a [351437] . In this model, the curvature 

# 3 It should be noted that we also need to take into account the normalization of the power spectrum of 
the curvature perturbations. 
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perturbation, up to the third order, is given by 



r Tc 

6 r 

r 

1296 



(T :k 



36- 



r 

72 

aacr 



6- 



+ (r 2 + 2r - 9) ^ 
18 (r 2 + 2r - 9) 



Sal 



(3r 4 + 10r 3 



r 3 

22r 2 - 54r + 135) -f 



(52) 



where r is the fraction of the energy density of the curvaton at the time of its decay and 
T a = dT /da and so on. In this model, the non-linearity parameters are written as 



/nl - - ( 3 



IT 



25 
54r 2 



36- 



r 2 i\ 



-2-r, 



18(r 2 + 2r-9)^ 



3r 4 + 10r 3 - 22r 2 - 54r + 135 



When r<l, /nl and #nl are related as 

2 
3 



r 2 r rr 

l 5 + 4 i i ^_ 18 1 



r-2 



rr 

-2^ + 3 



f 2 

J NL - 



(53) 



(54) 



(55) 



When the combination nVoyr 2 , is almost tuned to cancel the denominator in the coeffi- 
cient of /^L i n the right hand side, (/nl can be large even if / NL is 0(1) since the numerator 
is not necessarily canceled by such a choice of the functional form of V. 



3.2 Modulated curvaton 



Let us consider the case where the light field model contribution in Eq. (1511) comes from 
the modulated reheating mechanism. Assuming the quadratic potential of the light field 
and neglecting T aaa , in such case the total curvature perturbation is given by 

1 /_ . O _ Q\ i ( IVcr Tq- 

2 V ~ ~ f 2 

i_ 2r? 

36 



I — — I 5a* + 

1 3a, 6T ' 



+ 



81a 3 



-9r 2 + 



^ (3r-4r= 
+ 10r 4 + 3r 5 



3r a T a(J 



r a 



K 2 
5 at 



The non-linearity parameters are then given by 
/nl 



4r 
a. 



r 



r 



4r . 
+ ^(3-4r 



2r 2 



az 



(56) 



(57) 



50 / 4r 



fi'NL r — 

3 \ a* 1 



2Tj 

r 3 



+ + ^- (-18 + r + 20r 2 + 6r 3 ) 

9a 3 v 7 



(58) 
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From the above expressions, we find that for fine-tuned parameters we simultaneously 
have small /nl and large <7nl- F° r example, in case where T is given by Eq. (fT9|) and the 
parameters are set to be a = 1, r = A^ 1 and M/a, = A/4 with a big parameter A 3> 1, 
we can have a situation where 



/nl < 0(1), fi'NL = O(A), (59) 

by choosing (3 appropriately. Hence, we can simultaneously realize large #nl and /nl ^ 
0(1). 

3.3 Inhomogeneous end of hybrid inflation and curvaton 

In a similar way to the above discussion about the modulated curvaton mechanism, we 
consider the case where the fluctuations of the curvaton also induce inhomogeneous end of 
hybrid inflation discussed in 12.31 In this case, the total curvature perturbations are given 
by 



c= ( 2r | 1 pa, 



5a, 



+ 



+ 



9<x 2 
4 



(3r - 4r 2 - 2r 3 ) 



~ 9r + y 



2M 2 L a 
10r 4 + 3r 5 



4^2 



+ 



AM 



2r/cr 5-' 
1 f A ^(fal + 3v 



6M 2 t 



.9 



6^4 



5a 



+ 



1 palpal + Qy* 



Neglecting the Planck suppressed terms, the non-linearity parameters are given by 



fi 



NL 



5/2r | 1 pal 

3 V 3 ^ g 2 (f>cr 



50 /2r 



1 Pol 

Vcr g 2 4> 2 CT 

4 / 9 r 3 
— -9r 2 + — 
81 V 2 



- (3r - 4r 2 - 2r 3 ) + 
9 V ' 27/ 



1 Pa 2 Av 2 + Pa 2 " 



10r 4 + 3r 



<?V 4 r 



1 /V 4 (2/V 2 + 6u 



-,6^6 



(61) 



(62) 



In denominators in the above expressions, both terms are positive definite and hence in 
order to realize large onl the both terms must be much smaller than unity at least. For 
r < 1, in the numerator of the expression of /nl it is not possible to realize a cancellation 
to obtain the small /nl- Hence, the small /nl does not yield #nl large enough to be 
detected in this scenario. 
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3.4 Inhomogeneous end of thermal inflation and curvaton 

We can also consider the case where the light field contribution comes from the inhomo- 
geneous end of thermal inflation considered in 12.41 In such case, we have 



, ( 2r lg a \ 



( 3r _ 4r 2 _ 2r 3) + \( 9 - 

V ; 4 V 9 



9a 2 



9a 
2 



+ 



81o-3 



I + — 



) 9aa9a 



9' 
V 



Sal 
Sal 



The non-linearity parameters are then given by 



5/4r g 

Jnl — - 1 

3 \3a* 9 

50 / 4r g a 

S'nl = -zz 1 

27 \ 3a* p 



4r 
9oT 



( 3 _ 4r _2r 2 ) + (^-% 



2 \ i 



5 r 



!^ + 24 V ("18 + r + 20r 2 + 6r 3 ) 
# 2 2 7a 3 



(63) 
(64) 



(65) 

Similar to the case of the modulated curvaton case, for example, in case where 9 is given 
by Eq. ( 133]) and the parameters are set to be a = —1, r = A" 1 and M/a* = A/4 with a 
large parameter A, we have 

/nl < 0(1), Snl = O(A), (66) 

by choosing (3 appropriately. Hence, we can simultaneously realize large #nl and /nl < 
0(1). 

3.5 Modulated trapping and curvaton 

In case where a light field inducing modulated trapping mechanism behaves like curvaton 
in later dynamics, the curvature perturbation is given by 



_ ( 2r | 5e/T 
V3(x* 6a* 



(3r-4r 2 -2r 3 ) ' ^ 



9a 



81a- 3 



-9r' 



10r 4 + 3r 



8a 2 

5e/3 



Sal 



48a- 3 



Sat 



The non-linearity parameter are given by 



/nl = -(4r + 5e/3)" : 
25 

2nl = y(4r + 5e/3)" 
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4r (3 - 4r - 2r 2 ) + —eft 



16r 45 
— (-18 + r + 20r 2 + 6r 3 ) + — e/3 



(67) 
(68) 

(69) 
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In this scenario, due to the constraints on the parameters (3 and r as < (3 < 1 and 
< r < 1, it is hard to realize the large (?nl without violating the Planck constraint on 
/nl- 

4 Summary 

Following the Planck 2013 results [I], we discussed models of generating local-type non- 
Gaussianity where a light field other than the inflaton plays a main role of generating the 
curvature perturbations. First, we showed the constraint on model parameters for each 
light field model, by introducing the constraint on local-type /nl obtained in Planck 2013 
results: XXIV. By using the consistency relation between the non-linearity parameters 
/nl and <?nl for light field models as given in our previous paper [18], we also obtained the 
constraint on g^ and found that for simple light field models <?nl does not become large 
enough to be detected even in forthcoming Planck data and other future experiments, 
within the Planck constraint on /nl [3"BTPf2] . 

We also discussed the possibility of generating large #nl in light field models within the 
Planck constraint on /nl- We classified the possible models into two categories. One is to 
consider large third derivative of N, that is, N aaa in 5N formalism. Such examples include 
curvaton scenario with a self coupling and also modulated reheating (or inhomogeneous 
end of thermal inflation) scenario with non-trivial cubic dependence of a decay rate T (or 
a coupling g) on the modulus a which does not appear in the non-linearity parameter 
/nl- In such case, the consistency relation between / NL and #nl is no longer realized 
and hence we can realize the large (?nl within the Planck constraint. Another is to make 
the first derivative N a accidentally small, which leads to large #nl unless N aaa is also 
suppressed. Such accidentally small N a can be obtained in the case where the same light 
field u contributes to the curvature perturbations multiple times in different ways. As 
examples, we considered the cases where a modulus in a light field model also behaves like 
a curvaton at late epoch. We found that when we consider the modulated reheating or the 
inhomogeneous end of thermal inflation scenario as a light field model by taking fine-tuned 
parameters we can realize the large #nl due to the accidental cancellation. However, for 
the cases where the inhomogeneous end of hybrid inflation or the modulated trapping 
scenario is considered as a light field model we found it difficult to realize accidentally 
small N a and N aa and hence obtaining large <?nl within the Planck constraint is also 
difficult. Although information about the #nl would be a useful tool to discriminate the 
light field models in future observations, it seems to be difficult to realize the measurable 
g NL without fine-tuning. 
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